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We calculate the low temperature resistivity in low density 2D hole gases in GaAs heterostruc- 
tures by including screened charged impurity and phonon scattering in the theory. Our calculated 
resistance, which shows striking temperature dependent non-monotonicity arising from the compe- 
tition among screening, nondegeneracy, and phonon effects, is in excellent agreement with recent 
experimental data. 
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A number of recent density-dependent low tempera- 
ture transport measurements in dilute two dimensional 
(2D) n-Si MOSFET and p-GaAs heteroatructure systems 
have attracted a great deal of attentionEl because the ex- 
periments nominally exhibit a metal-insulator-transition 
(2D MIT) as a function of 2D carrier density (n). In ad- 
dition to this unexpected 2D MIT phenomenon (at this 
stage it is unclear whether the transition represents a 
true T = quantum phase transition (QPT) or a finite 
temperature crossover behavior) these measurements re- 
veal a number of intriguing transport propertieaj in di- 
lute 2D systems, such as a remarkable temperature de- 
pendence of the low density resistivity in the nominally 
metallic phase, which deserve serious theoretical atten- 
tion in their own rights irrespective of whether the 2D 
MIT phenomenon is a true QPT or not. 

In this paper we provide, a quantitative theory for 
one such recent experiment^ carried out in a low den- 
sity GaAs-based 2D hole gas. In our opinion, Ref. |^ 
represents a particularly important experiment in rela- 
tion to the 2D MIT phenomenon (although ironically 
no MIT is actually observed in Ref. ^ — even the low- 
est density data in Ref. ^ are entirely in the nominally 
metallic phase) because the ultra-pure samples used in 
Ref. ^ explore the 2D "metallic" regime of the highest 
mobility (i.e., the best quality or equivalently the low- 
est disorder), the lowest carrier density, and the low- 
est temperature so far studied in the context of the 2D 
MIT phenomenon. More spfijcifically, there have been 
suggestions and speculationall that the 2D MIT phe- 
nomenon is an interaction-driven QPT (the scaling the- 
ory of localizationB rules out a true localization transition 
in 2D disordered system) with the dimcnsionless Ts pa- 
rameter, which is the ratio of the interaction energy to 
the noninteracting kinetic energy of the 2D electron sys- 
tem, being the tuning parameter which drives the QPT. 
It is important to emphasize that increases as n de- 
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and therefore the 2D systems of 



Ref. y represent the highest (lowest) r^. [n) and conse- 
quently the most strongly interacting 2D systems exper- 
imentally studied so far in the context of the 2D MIT 
phenomenon. To be precise, values of the nominally 



"metallic" 2D hole regime explored in Ref. g go down 
to as low as = 26 (corresponding to the lowest hole 
density n — 3.8 x lO^cm^^ studied in Ref. ||) with no 
sign of an MIT whereas the other systema-studied in the 
literature exhibit the 2D MIT transitioro at critical r^ 
values as low as - 8 - 12 (Si MOSFETs) and 10 - 20 
(GaAs hole systems). The experimental results presented 
in Ref. |^ thus compellingly demonstrate that interaction 
(i.e., the parameter) is by no means the only (or per- 
haps even the dominant) variable controlling the physics 
of 2D MIT — disorder (and perhaps even temperature) 
also plays an important role. 

Our transport theory for the 2D hole system employs 
the finite temperature Boltzmann equation technique, 
which has earlier been successful in n-Si MOSFETscI 
and n-GaAs systemsdu. We include the following ef- 
fects in our calculation : (1) Subband confinement effects 
(i.e., we take into account the extent of the 2D system 
in the third dimension and do not assume it to be a 
zero-width 2D layer); (2) scattering by screened charged 
random impurity centers; (3) finite temperature and fi- 
nite wave vector screening through random phase ap- 
proximation (RPA) (actually we employ a slightly modi- 
fied versions of RPA, the so-called 2D Hubbard approx- 
imation, which approximately and rather crudely incor- 
porates the electron-electron interaction-induced vertex 
correction in the screening function which may be impor- 
tant at the low carrier densitieg3 being investigated — it 
turns out that our calculated resistance with the Hub- 
bard approximation is within 30% of the corresponding 
RPA results); (4) phonon scattering!. The effects we ne- 
glect in our theory are (1) all localization and multiple 
scattering corrections; (2) inelastic electron interaction 
effects — in fact, all effects of electron-electron interac- 
tion are neglected in our theory except for the long range 
screening through RPA and (approximate) short-range 
vertex correction through Hubbard approximation. 

Our calculations are similar to the onesEl we recently 
carried out for electron inversion layers in n-Si MOSFETs 
with two important differences; (1) we include the full 
hole density in the current calculations without subtract- 
ing out any critical density as done in Ref. H — this is, in 
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fact, consistent with our Si MOSFET calculations since 
the critical density in Ref. || must be extremely small, 
and in any case SdH measurements carried out in Ref. 
[2] show that all the carriers are "free" and participating 
in the conduction process; (2) we include phonon scat- 
tering effects in the current calculations because phonon 
scattering is significant for GaAs holes already in the 
T = 1 — lOK temperature range whereas phonon scatter- 
ing is neghgibly small in n-Si MOSFETs in the 1 — lOK 
temperature range. Details of phonon scattering calcula- 
tions are given in Ref. ^ — the essential point is that the 
phonon resistivity is proportional to T for T > IK and is 
negligibly small in the low temperature Bloch-Griineisen 
regime. 

Our calculated resistivity for 2D holes in GaAs struc- 
tures is shown in Figs. 1 and 2 for two different types of 
2D quantum confinement: Square well (Fig. 1) and het- 
erojunction (inversion layer type approximately "trian- 
gular") confinement (Fig. 2). The qualitative results for 
the two kinds of confinement are, as expected, very sim- 
ilar (although the actual quantitative resistance values 
depend on the nature of confinement since the scattering 
and screening matrix elements are strongly confinement 
dependentu through the wavefunction spread normal to 
the 2D confinement plane). The resistivity can be written 
as p{T) = pq + PimpiT) + Pph{T), where po = p{T -> 0) 
is the residual resistivity arising entirely from (screened) 
charged impurity scattering in our theory (for a weakly 
localized system po diverges logarithmically as T — > 0, 
our theory is valid above the crossover temperature scale 
for weak localization to set in — no indication for the 
expected In T weak localization divergence is seen in the 
experimental data of Ref. [2] down to the lowest reported 
measurement temperature, SSmif). /5p?i(K) is the resis- 
tivity contribution by phonon scatterings which could 
be quite significant for 2D holes in GaAs already in the 
1 — IQK temperature range. Finally, Pimp{T) is the tem- 
perature dependent part of the charged impurity (i.e., ran- 
dom disorder) scattering contribution to the resistivity, 
i.e., po + Pimp(T) = Pi is the total impurity contribution 
to the resistivity. We note that po, which sets the overall 
resistivity scale [by definition, both pimp{T) and pph{T) 
vanish as T — > 0] in the problem, is determined by the 
amount of the random disorder in the system which is in 
general unknown. The amount of random disorder (and 
consequently po) depends on the strength and the spa- 
tial distribution of all the impurity scattering centers in 
the system. We parameterize the charged impurity den- 
sity, assuming them to be randomly distributed static 
Coulomb charged centers interacting with the 2D carri- 
ers via the screened Coulomb interaction. We adjust the 
charged impurity density (assumed to be randomly dis- 
tributed in our calculations) to get agreement between 
theory and the experimental data — thus the scale po is 
essentially an adjustable parameter in our theory since 
the actual impurity distribution in the 2D systems of in- 
terest is simply not known. We emphasize, however, that 
the charged impurity density needed in our theory to 
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FIG. 1. The calculated 2D hole resistivity p(T) for sym- 
metric square well systems corresponding to the hole densi- 
ties n = 0.65, 1.07, 1.63, 3.26, 4.80x lO^^cm"^ (from top to 
bottom) with random impurity densities Ni = 0.7, 0.75, 0.8, 
2.0, 3.5x lO^^cm"^, respectively. In inset the sample configu- 
ration is shown schematically. In this calculation we use the 
parameters, ~ 300A, di = 200^4, and d2 = 50A. We use 
a very small random impurity density, N^, = 3 x lO^^cm"^, 
in the GaAs layer itself which is consistent with the extreme 
high quality of the samples in Ref. [2]. Some representative 
experimental data points from Ref. [2] are shown (the actual 
random disorder in the experimental samples is unknown). 

obtain agreement between our calculations and the ex- 
perimental data for po are reasonable. 

Before discussing our results we make three salient re- 
marks about our calculation and model. First, we ne- 
glect scattering by interface roughness, alloy disorder, 
etc. in our calculation (including only charged impurity 
scattering in the theory) since it is well-known that the 
dominant low temperature resistive mechanism in high 
quality GaAs structures arises essentially from charged 
impurity scattering (it is straightforward to include addi- 
tional scattering mechanisms in our calculations with the 
unpleasant complication of having additional unknown 
parameters, such as the interface roughness strength, in 
the theory — our choice is to keep the number of un- 
known adjustable parameters at a minimum by assum- 
ing that all of the random disorder scattering is caused by 
randomly distributed charged impurity scattering which 
should be an excellent approximation for the extreme 
high quality GaAs samples used in Ref. [2]). Second, the 
Matthiessen's rule, which is implicitly assumed in sepa- 
rating out Pi(T) and pph{Tl, is known to be not strictly 
valid at finite temperaturesu because different scattering 
rates do not simply add in the total resistivity. It is im- 
portant to emphasize, however, that we do not assume 
the Matthiessen's rule in our theoretical calculations, and 
Eq. (1) is written down simply as a rough guide for 
qualitative discussion. In any case, the deviation from 
Matthiessen's rule is of the order of 30% or less, which is 



2 



not of much consequence for our discussion. Finally, the 
third remark we make is regarding our use of the single 
scattering Born approximation in our Boltzmann theory 
(neglecting all multiple scattering effects), which can be 
justified by noting that our calculated resistivity (and 
the corresponding experimental resistivity measured in 
Ref. [2]) always satisfies the weak scattering condition of 
kpl ^1 — in fact, our results are restricted to kpl > 3 
even in the worst situation (for our highest resistance re- 
sults) . We therefore believe that the Born approximation 
may not be a poor approximation for our problem. 

In Fig. 1 we show our calculated 2D hole resistivity 
for symmetric square well systems corresponding to the 
sample of Ref. [2]. The actual sample configuration is 
shown schematically as an inset in. Fig. 1. We also show 
some representative experimentala results (from Fig. 2 
of ref. [2]). We emphasize that the quantitative agree- 
ment with the data of ref. [2], while being certainly in- 
dicative of the essential validity of our theoretical ap- 
proach, should not be taken too seriously — it is cer- 
tainly not the feature of our theory we would focus on, 
particularly since the random impurity distribution in 
the experimental samples is unknown. It is the overall 
striking qualitative similarity between our microscopic 
theory and the experimental dataa which deserves at- 
tention. This is particularly so because the density and 
temperature dependence of the measured resistance in 
ref. [2] shows a throughly nontrivial non-monotonic be- 
havior which is completely reproduced in our calcula- 
tions. This striking non-monotonicity in p{T), at lower 
carrier densities, arises from a competition among three 
mechanisms: Screening, which is particularly important 
at lower T; nondegeneracy and the associated quantum- 
classical crossover for T >Tf {= Ep/kB, the Fermi tem- 
perature) which was discussed in ref. ^ in the context 
of n-Si MOSFETs; and phonon scattering effect which 
is negligible below IK, but starts becoming quantita- 
tively increasingly important for T > IK. The Fermi 
temperature for the 2D hole system can be expressed as 
Tp = 0.64(n/10^°)X where n is the 2D hole density mea- 
sured in units of 10^°cm~^. Thus for n = 4.8 x lO^^cm"^ 
between n — 0.65 x lO^^^cm"^ in Fig. 1 Tp varies be- 
tween 3K and OAK. This makes the quantum-classical 
crossover physics particularly significant for the results 
of Ref. [2] as was already noted by the authors in Ref. 
[2]. 

At higher densities (the bottom two curves in Fig. 1) 
the quantum-classical crossover effects are not particu- 
larly important because phonon scattering becomes im- 
portant before the classical behavior! p can show 
up, and the system makes a transition from the quantum 
regime to the phonon scattering dominated regime — the 
fast rise in p{T) at high T in Fig. 1 is the phonon scat- 
tering effect. At low enough densities, however, phonon 
scattering effects are absent (because phonons are frozen 
out in the low temperature Bloch-Griineisen rangeB) at 
the quantum- classical crossover point which occurs at 
very low temperatures around T < Tp < IK (the top 
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FIG. 2. The calculated resistivity for a heterostructure in- 
version-layer-type "triangular" confinement 2D (a) hole gas 
and (b) electron gas for carrier densities n — 0.38, 0.65, 
1.07, 1.63, 1.93, 3.26, 4.15, 4.80, 8.66 xlO^^cm"^ (from 
top to bottom) with random impurity densities Ni = 0.3, 
0.3, 0.3, 0.3, 0.3, 0.9, 1.1, 1.1, 5.0 xlO^^cm-^, respectively. 
The inset in (a) show the sample configuration. We use 
di — 250A and d2 = lOOA. The impurity density in GaAs is 
7V„ = 3 X lO^^cm"^ 

two curves in Fig. 1). In these low density results one 
can see p{T) increasing pwith T at lower temperatures 
due to screening effectsEl, then the quantum-classical 
crossover occurs at the intermediate temperature regime 
around where nondegeneracy effects make resistivity 
decreasda as p '~ T''-^; eventually at higher temperatures 
(T > IK) phonon scattering takes over and p{T) in- 
creases with T again. At higher densities Tp is pushed 
up to the phonon scattering regime, and the quantum- 
classical crossover physics is pre-empted by phonons so 
that non-monotonicity effects are not manifest. 

The non-monotonic behavior of p{T) as a function of 
n and T is made more explicit in Fig. 2(a) where we 
show our calculated resistivity for the same density and 
temperature range as in Fig. 1 for a heterostructure 
inversion-layer-type "triangular" confinement 2D hole 
gas, separating out the pure impurity scattering contri- 
bution (i.e., the dashed curves in Fig. 2(a) leave out 
the phonon scattering contribution completely). First, 
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we note that the resistivity results in Fig. 2(a) are very 
similar to those in Fig. 1, indicating that the transport 
behavior seen in ref. [2] is the generic behavior of a low 
density 2D GaAs hole system, and does not arise from 
any particular feature of the square well samples used in 
ref. [2]. Second, the interplay of screening (low temper- 
ature), phonons (high temperature), and nondegeneracy 
(high temperature and low density) is manifestly obvious 
in Fig. 2(a): the intriguing low density non-monotonicity 
in the observed p{T) clearly arises from the fact that 
both screening and phonon scattering mechanisms give 
rise to a p(T) monotonically increasing with T (at low 
temperature for screening, and at high temperatures for 
phonons), but nondegeneracy effects produce a p{T) de- 
creasing with T for T > Tp- Since phonon scattering is 
the dominant temperature dependent scattering mecha- 
nism in GaAs holes for T > IK, the non-monotonicity 
can show up in any significant way only if Tp < IK, 
which is precisely the experimental observation. 

As an interesting comparison we show in Fig. 2(b) the 
calculated p{T), without any phonon scattering, for the 
same densities (and impurity scattering parameters) as 
in Fig. 2(a) for a 2D electron inversion layers confined 
in a GaAs heterostructure (i.e., the only difference be- 
tween the results for Fig. 2(a) and Fig. 2(b) is that the 
GaAs electron mass has been used in the calculations 
corresponding to Fig. 2(b) rather that the hole mass. 
The neglect of phonon sfiattering is justified by the fact 
that phonons contributeQ significantly to GaAs 2D elec^ 
tron resistivity only for T > lOK — in fact, inclusion^ 
of appropriate phonon scattering would produce results 
indistinguishable from the results shown in Fig. 2(b) 
(i.e. upto 5K). The difference between the results of 
Figs. 2(a) (holes) and 2(b) (electrons) is striking: there 
is essentially no observable (on log scale) temperature 
dependence at low temperatures in the 2D electron re- 
sistivity in GaAs heterostructure down to 2D densities 
as low as n = 0.38 x 10^"cm~^. This essential temper- 
ature independence of low temperature electronic resis- 
tance in highL quality GaAs heterostructures, which is 
a well-knowrjj experimental fact, arises from the weak 
screening property (associated with its low effective mass 
and the associated small electronic density of states) of 
2D electrons in GaAs heterostructures compared with 
higher mass 2D holes in GaAS or 2D electrons in Si MOS- 
FETs. This weak screening behavior of GaAs electrons 
precludes any strong temperature dependent p{T) even 
at very low carrier densities (and temperatures). The 
quantum-classical crossover phenomenon, however, still 
occurs around T ~ Tp, leading to a p{T) ^ T^^ for 
T > Tp, which is manifestly obvious in Fig. 2(b), par- 
ticularly for lower densities. Note that the Fermi tem- 
perature in Fig. 2(b) corresponds to Tp — 4:.l{n/lQi^^)K 
with n being the 2D electron density in Fig. 2(b) mea- 
sured in units of lO'^^cm"^. Thus the Fermi temperature 
in Fig. 2(b) ranges from l.bK (top curve) to ib.bK (bot- 
tom curve). We note that the decreasing p{T) at higher 
T in Fig. 2(b) arises not only from a quantum to classical 



crossover (which is the dominant effect at lower densities 
when Tp is low), but also from the finite temperature 
Fermi surface averaging in a degenerate quantum sys- 
tem. It is easy to show that the Fermi surface averaging 
effect at finite temperatures, by itself, always leads to 
a finite temperature resistivity which decreases weakly 
with temperature (even in the T — > limit) - in fact, 
this effect by itseff leads to p{T) « po[l - 0{T/Tp)% 
and can only be observed if the temperature dependent 
screening effects are unimportant. This effect was first 
observed in 2D electro»s in GaAS heterostructures more 
than fifteen years agoM. 

To conclude, we have developed a theory for the low 
temperature transport properties of 2D holes and elec- 
trons confined in low density and high mobility GaAs het- 
erostructures. Our theory includes temperature depen- 
dent screening of impurity scattering and phonon scat- 
tering effects. Agreement between our theory and ex- 
periment suggests that screening and impurity scatter- 
ing effects play an essential role in determining much of 
the intriguing temperature and density dependent trans- 
port properties in 2D systems, and that random disorder 
(mostly arising from charged impurity scattering) is an 
important ingredient in the physics of low density 2D 
systems. 
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